SYMMETRY OF BIPOLARON BOUND STATES 
FOR SMALL COULOMB REPULSION 
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Abstract. We consider the bipolaron in the Pekar-Tomasevich approximation and 
address the question whether the ground state is spherically symmetric or not. Nu- 
merical analysis has, so far, not completely settled the question. Our contribution 
is to prove rigorously that the ground state remains spherical for small values of the 
^1 ' electron-electron Coulomb repulsion. 

S3 ' 



1. Introduction 



On 

In this paper we shall be concerned with properties of the bound state of two 
Ph. polarons. We do this in the context of the Pekar-Tomasevich model (T5J [T6] of the large 

polaron, which, in turn, is based on Frohlich's polaron model [5]. In the latter model 
two electrons interact with a quantized electric field generated by the displacement of 
the nuclei in a polar lattice. There are two coupling constants in Frohlich's model. 
The coupling to the field, a, and the Coulomb repulsion among the electrons, denoted 
by U. Pekar's approximation is to assert that the wave function is a product of a 
two-particle electron wave function %p times a field function $. After eliminating the 
field one is led to Pekar's energy expression for ip, 

£u[1>]= ft (\V x ^\ 2 + \V y i;\ 2 + —^\^Adxdy-2aD\p il ,p^. (1.1) 

J Jr 3 xi 3 V \ x ~ y\ J 

The electron coordinates are x and y and the electron spin does not appear explicitly, 
except that ip is symmetric for the ground state, which is a singlet state. (The reason 
that it is a singlet is that the ground state is a positive function and must, therefore, 
be symmetric - an observation that goes back to Wigner many years ago.) 
In ( 11. II) p^ denotes the electron density, given by 

Pi>{x) = / \ip(x,y)\ 2 dy + \xp(y, x)\ 2 dy , 

JM.3 J R 3 

and D[p, p] is the Coulomb energy of a charge distribution p, 

D\P,P} = 1 II P \ X)p{X j} dxdx>. 



2_ 

Note the minus sign in (jl.ip ; the induced interaction is attractive. 



© 2012 by the authors. This paper may be reproduced, in its entirety, for non-commercial 
purposes. 



2 



R. L. FRANK, E. H. LIEB, AND R. SEIRINGER 



The Pekar energy is 

e v = MiSuM : ^ e H\R 6 ) , \\i>\\ = 1} , (1.2) 

with H^ll denoting the L 2 norm of if) and H 1 denoting the Sobolev space, i.e., square- 
integrable functions whose gradient is also square-integrable. This energy is more than 
just an approximation, for it is asymptotically exact as a and U tend to infinity with 
U/a fixed. This is stated in [13] . following the technique of [IB]; see also [2]. 

A benchmark for the bipolaron problem is the energy of a single polaron. It is 
defined in the manner of equation (11.11) 

m = [ W\ 2 dx-2aD[\ij\ 2 ,\ij\ 2 }, 

and 

e = M{£&] : V e ^(R 3 ) , IMI = 1} • (1-3) 
It is known that there is a minimizing ip for this single polaron problem and that it is 
unique, up to translations in R 3 and multiplication by a constant phase [10]. Since it 
is unique, it is a radial function. 

By scaling we can always reduce to the case a = 1/2 and we shall do so henceforth. 

There is a considerable literature on the subject of rotation invariance of the bipo- 
laron energy minimizer, usually formulated in the language of 'one-center bipolaron 
versus two-center bipolaron'. The analyses are all based on variational calculations. 
While there seems to be general agreement that the one-center bipolaron has the lower 
energy, it is not completely clear that a more sophisticated variational treatment will 
preserve rotational symmetry, especially near the value of U where the bipolaron ceases 
to be bound. We have shown rigorously that there is such a critical U c [3] . Numerical 
variational data seems to indicate that the critical U c is rather close to 1, namely, 
U c ~ 1.15. This tells us that binding is a delicate matter and, indeed, the existence 
of such a U c was an open question for some time. 

A minimizer also exists for a bipolaron provided the energy is below the energy of 
twice the single polaron energy [8]. Interestingly, the bipolaron has a minimizer with 
finite radius at the critical value U = U c [3J. (The same holds for a helium atom for 
the critical value of the nuclear charge.) The existence of a minimizer implies that the 
translation invariance of the Pekar minimization problem is broken. It is, therefore, 
not out of the question that the minimizer for two polarons might break rotational 
symmetry as well in order to lessen the Coulomb repulsion. 

The value of U determined by physical electrostatic considerations is always U > 1. 
Nevertheless, one can consider the mathematical question for small, but positive U 
and ask whether there is a possible lack of rotational invariance in that case. After 
all, a rotating object like the earth becomes oblate even for the smallest amount of 
rotation. 

In this paper we will prove that there is no breaking of rotational symmetry for small 
U. Our strategy for proving the lack of symmetry-breaking for small U is based on 
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the following consideration. For any small U the minimization problem, restricted to 
rotation invariant functions, has an energy minimizer, as we shall prove. The question 
is whether there is better minimizer without rotation invariance. Since U is small, 
both minimizers would have to be very close to the unique (up to translations) U = 
minimizer, and so we can discuss the existence of a symmetry breaking minimizer 
by means of rigorously controlled perturbation theory. It is evident that the benefit 
of symmetry breaking to the repulsive energy will be proportional to — 5 2 , where 5 
measures the non-sphericity, but this contribution is multiplied by U. On the other 
hand, the increase in the rest of the energy is presumably also of the form AS 2 , where 
A > is some U independent number. If A > then 5 wants to be zero for small U. 

The problem with this argument is that A could be zero, in which case the —U5 2 
energy would always win, no matter how small U is. Most of what we do in the paper, 
from the mathematical point of view, is to show rigorously that A is not zero, and 
thus there is no distortion for small U. 

Even utilizing the result that A > for the single polaron [7] , the proof given here 
for the bipolaron will not be a short one. One of the complexities faced in this proof 
is the fact that a simple translation is a distortion that costs no energy, i.e., there are 
zero modes. While these are physically trivial distortions it is not a trivial matter to 
separate their contribution, mathematically, from the relevant ones. 

Theorem 1. There is a U s > such that for all U < U s the minimizer of Sjj is 
unique up to translations and multiplication by a constant phase. In particular, after 
a translation it is rotation invariant, that is, ?J)(lZx,lZy) = ip(x,y) for any x,y EM. 3 
and any 71 G 0(3). 

It remains an open problem to decide whether the ground state ceases to be rotation 
invariant for U close to the critical value U c . 

The rest of this paper has two parts. In Part A, we reduce the proof of Theorem[TJto 
a problem in second-order perturbation theory, namely the question of the distortion 
coefficient A. In Part B, we complete the proof of Theorem [1] by showing that A is 
positive and that zero modes play no important role. 

Acknowledgments. We are grateful to Herbert Spohn for making us aware of this 
problem. Partial financial support from the U.S. National Science Foundation through 
grants PHY-1068285 (R.F.), PHY-0965859 (E.L.) and the NSERC (R.S.) is acknowl- 
edged. 



2.1. Some preparations. Step 1. In searching for the minimum in fll.2p we can 

confine our attention to non-negative, symmetric (i.e., ip(x,y) = tjj(y,x)) functions. 
The reason is that we can replace any ip by 



2. Proof of Theorem [Q Part A 




(2.1) 
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The potential energy terms remain the same and the kinetic energy term does not 
increase [HJ Thm. 7.8]. 

Step 2. An important step is to reformulate the energy minimization problem in 
the following way: Define 

£u[il>,$]:= II (\V x ^\ 2 + \V y ^\ 2 ~^xM\ 2 -^yM\ 2 + -^-\i;\ 2 ) dxdy 

j jk 3 xr 3 v \ x ~y\ j 

+ — I \V$\ 2 dx If \^\ 2 dxdy (2.2) 
for if) G H 1 ^) and $ G if 1 (R 3 ). Then 

^[^$]>ii^n 2 ^[ii^irv] (2.3) 

with equality if and only if $ = * Pv/lhAlh 

One advantages of this reformulation is that one can see immediately that there is a 
unique (up to translations) minimizer for the U = problem. For a given $ we have, 
in this Schrodinger minimization problem for two independent particles, so the 

best if) is a product, ^(x,y) = f(x)f(y), in which case p$ = 2\f\ 2 and the optimal $ 
is 2|x| _1 * \f\ 2 . The problem therefore becomes 

e = inf {2 £ \Wf\ 2 dx - 4D[|/| 2 , |/| 2 ] : ||/|| = ll , 

which coincides, up to a rescaling, with problem ( 11. 3D . We conclude that eo = 8e and 
that the minimizing / is unique up to translations and multiplication by a constant 
phase [10] . 

Another advantage of this reformulation is that one sees that for any U > an 
optimizer ip, if it exists, is the ground state of a two-body Schrodinger operator. 
Therefore, ip is a multiple of a strictly positive function [TT| Thm. 9.10]. This, in turn, 
implies that i/j is symmetric (i.e., ip(x, y) = ip(y, x)), because otherwise replacing ip by 
(12.11) would strictly lower the kinetic energy [TTJ Thm. 7.8]. 

Step 3. In addition to the global minimization problem (11.21) one can define the 
rotationally symmetric minimization problem, that is, 

^symm _ /^jy/j • ^ (z Jj 1 (M?) , ||^|| = 1, if) rotation invariant with respect to 0} . 

(2-4) 

N.B.: From now on 'symmetry' refers to 'rotation symmetry' and not to symmetry in 
x and y. We recall that the rotation invariance of ip means that ip(TZx, TZy) = ip{x, y) 
for any x,y G M 3 and any 71 G 0(3). The density of such ip is, of course, radial 
and, by Newton's theorem, its potential p^ * is a symmetric decreasing function 
bounded by 2|x|~ 1 . 

Of course, eu < e s J mm , and our goal in this paper is to investigate whether equality 
holds. We collect some properties of these energies. Both ejj and e^f mm are non- 
decreasing, concave functions of U (as infima of non-decreasing, linear functions). 
Because of Step 2, e = e s ymm for U = 0. Moreover, simple trial function arguments 
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show that eu < 2e and e s J mm < e for all U . Lewin [8] has shown that the infimum eu is 
attained provided eu < 2e. In the appendix of this paper we shall prove an analogous 
result for the rotation invariant problem, with a different condition, however, namely, 

symm 

e u < e - 

Proposition 2. If e^ mm < e, then the infimum in (12.41) is attained. 

The reason for the discrepancy between 2e for eu and e for e s J mm [ s that spherical 
symmetry prevents the formation of two more or less separate polarons. In other 
words, the second polaron density has to be far away from the first, and in a radial 
shell, which makes it impossible to retain an energy 2e with an essentially unbound 
pair of polarons, both of which are spherically symmetric with respect to a common 
center. 

Proposition [2] shows the fact that for some values of U the rotation invariant mini- 
mizer is not the true minimizer. A rotation invariant minimizer is necessarily a critical 
point of the Pekar-Tomasevich functional. If it is a 'false' minimum (i.e., its energy 
is bigger than the true ground state energy), its existence can possibly lead to com- 
putational difficulties for the true minimizer. 

In the appendix we show two things about the rotation invariant minimization 
problem: (1) There is a critical constant U^ ymm with 1 < U^ ymm < 4 such that 

e symm < g fo r £/ < Jjsymm &nd ^jmm = g fo r £/ > Jjsymm. ^ p roposition gj 

(2) The rotation invariant minimizer, if it exists, is a function of \x\, \y\ and t = 
x ■ y/\x\\y\. We will show that for fixed |x| and \y\ the minimizer is non- increasing as 
a function of t. That is, the two particles try to avoid each other; see Proposition |9j 



2.2. Beginning of the proof of Theorem [TJ After these preparations we are now 
ready to give the first part of the proof of Theorem [TJ Our strategy is as follows: 
Let U n > be a sequence such that U n — > 0. For every sufficiently large n there is a 
global minimizer ip n and a rotation invariant minimizer (p n corresponding to fll.2[) and 
(12.41) with U = U n . We shall prove that for all large n, cp n and tp n coincide up to a 
translation and a constant phase. 

This clearly implies the theorem. Indeed, if the theorem were not true we could find 
a sequence U n tending to zero, and associated minimizers ip n and ip n which are not 
translates or multiples of each other. By what we are going to show, however, they 
are both translates and multiples of a rotation invariant minimizer cp n if n is large, 
which is a contradiction. 

Thus, from now on we fix a sequence U n > tending to zero. We will only consider 
n sufficiently large such that there are a global minimizer ip n and a rotation invariant 
minimizer ip n corresponding to (11.21) and (12. 4p with U = U n . By Steps 1 and 2 
above we may assume that ip n and tp n are positive and permutation symmetric (i.e., 
ip n (x,y) = i[) n (y,x) and similarly for ip n ). Since the global minimization problem is 
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translation invariant, we may translate ip n in such a way that 



«n := // ip n (x,y)ip n (x } y)dxdy = max / / ijj n (x - a,y - a)ip n (x,y) dx dy . 
7jr3 xK 3 am 3 j7r3 X]R 3 

(2.5) 

With this normalization, our goal is to show that ip n = ip n for all large n. 
We decompose 

where, according to the definition of a n and to the maximizing property in (12. 5p . 

jnfn dx dy = and / / j n e • (V x + V y ) i Pn dxdy = for any e G S 2 . 

1 J JR3xR3 

(2.6) 

The derivative term comes from differentiating the last double integral in (I2.5P with 
respect to a. 

We next claim that 

Ibnllif 1 — ► as n — )• oo . (2.7) 

Indeed, both ip n and <p n are minimizing sequences for the eo problem. Since the 
minimizer for this problem is unique up to translations (see Step 2 above), the results 
of [8] imply that 

ipn{- ~ an, • ~ a n ) / ® / and <p n {- - b n , • - b n ) -> f ® f in if^R 6 ) 

for some sequences a n and 6 n in M 3 . Here / ® / is the eo minimizer with / chosen 
spherically symmetric about the same origin that we fixed to formulate the problem 
(12. 4ft . Since the density of <p n is spherically symmetric, one easily concludes that 
b n — > and, therefore, <p n — > f ® / in if 1 as n — )■ oo. This, together with the 
maximizing property of ot n in ( 12. 5ft . implies also that a n — )■ and — )■ / <g> / in if 1 
as n — i- oo. Thus a n — >• 1 and j n — — ctnfn — > in if 1 , as claimed in ( 12. 7ft . 

We now expand the energy of ip n to second order in Hjnlltf 1 - A simple but tedious 
computation shows that 

SuM = £uM + (jn: L n Jn) + 0(\\j n \\ 3 Hl ) (2.8) 

with the linear operator 

-A s - A y + Un - P(pn * |^ - Pipn * J- - ^ J - 4A^ n . 
p — 2/| fI |2/| / 

Here p ¥ , n * | xz? | 1 is an abbreviation for (p Vn * \ ■ | _1 ) (x), we introduced fi n = £u„[(Pn] — 
D[ptp n , pp n ] and X^ is the integral operator on L 2 (M 6 ) with the integral kernel 

Y ( / i\ _ ( Pn(x : y)ip n (x',y') t (p n (x,y)ip n (x',y') 



\y-y 
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In order to see that there is no linear term in j n in the computation, we used the 
Euler-Lagrange equation for (p n , that is, 

\-A x - A y + Un - p Vn * - p Vn * J ip n = p n ip n . (2.9) 
V \x — y\ \x\ \y\J 

(From the minimizing property of (p n one obtains this equation only when integrated 
against rotation invariant functions but, since both sides are rotation invariant func- 
tions, it is true even when integrated against any function.) We also used the fact 
that 

«n=(l-||jn|| 2 ) 1/2 = l-^||jn|| 2 + 0(||jn|| 4 ). 

In deriving (12. 8 j) we applied standard estimates to bound all terms of higher than 
quadratic order by a constant times Hjnlljyi- 

The way forward is now clear: The quantity A referred to in the introduction can 
be identified as (j, L n j)/\\j\\ 2 for j's satisfying (12. 6p . and we need to prove that this 
is non-zero. We will prove that there are constants N > 1 and c > such that 
(j,L n j) > c\\j\\ 2 for all n > N and all j satisfying the orthogonality conditions (12. 6ft . 
Then (ESJ) will imply that 

SuM > S Un [p n ] + c\\ 3n \\ 2 + 0(\\j n \\ 3 m ) . 

Since ip n is, by assumption, an energy minimizer, we have necessarily j n = for all 
large n, which means ipn = fn, as we intended to prove. 

Thus we are left with proving a lower bound on (j, L n j). We show this perturbatively 
by analyzing the U = case. 

2.3. The Hessian. Instead of working directly with the operator L n , it is more con- 
venient to work with a closely related operator, namely the Hessian H n of E\j n at tp n . 
That is, for any normalized, real- valued j G if 1 (R 6 ) we define 



de 2 



e=0 



¥n + ej 



(j,H n j). (2.10) 



y/l + 2e{<p n ,j)+e 2 
A similar computation as before shows that 

H n = L n + \k n )(ip n \ + \(p n )(k n \ + PnWn) (fnl , 

where 
and 

kn = -2 ~ Ay + \ x U -y\ ) Vn ~ 2Pv>n * l x l~ lv9n ~ 2pv>n * \y\~ 1( Pn 

The expressions for (3 n and k n can be somewhat simplified using equation (12. 9p for ip r 
but we will not need this. The only thing that is relevant for us is that (j, H n j) = 
(j,L n j) if j satisfies the first orthogonality condition (I2.6p . 
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We collect two facts about the operators H n . First, H n (which commutes with 
angular momentum since ip n is rotation invariant) is non-negative on the subspace of 
angular momentum zero. This follows from the minimizing property of ip n . Second, 
the functions <p n and e^Vas+Vj,)^, e G § 2 , are in the kernel of H n . For ip n , this follows 
immediately from the Euler-Lagrange equation ( 12. 91) for ip n , and for e • (V x . + V y )ip n) 
this follows by differentiating the equation (I2.9p for ip n (x + te,y + te) (which is the 
same as that for ip n ) with respect to t at t — 0. 

The following proposition says that, for small U n , H n is strictly positive away from 
the zero modes found above. 

Proposition 3. There is a number c > and an N > 1 such that for all n > N 
and for all j which are orthogonal to (p n and to e • (V x + V y )ip n , e G S 2 , the following 
estimate holds: 

(j,H nJ )>c\\j\\ 2 . 

Equivalently, if P n denotes the projection onto the four- dimensional space spanned by 
(p n and e ■ (V x + V y )(p n , e G S 2 , then for n > N 

P L H P _L > cP L 

Note that, as remarked above, (j,H n j) = (j,L n j) if ((f n ,j) = and, therefore, 
Proposition |3] concludes the proof of Theorem [TJ □ 

Since f n — > f ® f in if 1 (IR 6 ) as n — > oo, the operators H n tend to the corresponding 
operator at U = in norm resolvent sense. In particular, the eigenvalues converge 
and, therefore, it suffices to prove Proposition [3] in the case {7 = 0. This is the topic 
of the next section. 

The reason that we have chosen to work with the operator H n instead of L n is that 
the operator L n has a negative eigenvalue. (This follows from the variational principle 
since ((p n ,L n ip n ) = -4((p n X Vn ,(p n ) < by §2Jfl).) The positivity of (j,L n j) asserted 
in Proposition |3] therefore crucially relies on the orthogonality condition (12. 6p . This 
condition is not easy to use, however, since ip n is not an eigenfunction of L n . In 
contrast, it is an eigenfunction of H n . We also note that H n is the operator that 
automatically takes care of the normalization condition, without any reference to 
orthogonality, and, therefore, is most directly connected to the coefficient A mentioned 
in the introduction. 



3. Proof of Theorem [H Part B 

As explained in Step 2 of the previous section, at U = we have <p(x, y) = f(x)f(y), 
where / is a radial decreasing function on M 3 with j f 2 dx = l and 

hf = 0, h = -A-2f*±--^. (3.1) 

\X\ — 
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Here, fx = e -AD[f 2 , f 2 } = 8e-AD[f 2 , f 2 ]. The Hessian of £ at /<8>/ is the operator 
on L 2 (R 6 ) given by 

H = L+ \r ® / + / ® r)(f ® /| + |/ ® /)<r ® / + / ® r| + ® /)(/ ® /| , 

where 

L = h x + h y — 4X/0j 
and where /3 is a constant and r is a function in L 2 (R 3 ). More precisely, 

(3 = 8{(f,-Af)-6D[f 2 J 2 ]) 

and 

r = -2 (-A -4/ 2 * (x^ 1 ) /. 

The main result of this section is the following U = analogue of Proposition [3l We 
note that the three-dimensional space spanned by e • (V x + V y )<£>, e G S 2 , coincides 
with the space spanned by / ® f'Y 1>m + f'Yi im ® /, m = —1,0,1, where Fi >TO are 
spherical harmonics of degree one. 

Proposition 4. There is a number C > suc/i that for all j G i/ 1 (IR 6 ) with j(x,y) = 
j(y, x) that satisfy 

(/ ® /, J) = (/ ® + ® /, j) = /or all m = -1, 0, 1, 

i/ie following estimate holds: 

(j,Hj)>C\\j\\ 2 . 

Equivalently, if P denotes the projection onto the four- dimensional space spanned by 
f®fandf® f'Y 1>m + f'Y 1>m ®f,m = -1, 0, 1, then 

P ± HP ± > CP ± . 

We shall deduce this from a result about the Hessian of the one-particle functional, 
which we discuss in the following subsection. 

3.1. The one-particle Hessian. We recall that the function / minimizes the one- 
polaron functional 

/ |W| 2 dx-2D[H 2 ,H 2 ] (3.2) 

and that the corresponding Euler-Lagrange equation reads hf = with h from ( 13. lft . 
Moreover, the Hessian of the above one-polaron functional at / reads 

^ 1 ) = L( 1 ) + |r)(/| + |/)(r|+ 7 |/)(/| 

with 

= h-Ax } . 

Here xj is the operator on L 2 (R 3 ) with the integral kernel f{x)\x — x'\~ l f{x') and r is 
the same function as in the expression of the two-particle operator H. The following 
theorem is the one-particle analogue of Proposition HI Its proof relies heavily on 
previous work of Lenzmann [7] . 
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Proposition 5. There is a number C > such that for all j G i/ 1 (IR 3 ) that satisfy 

(f,j) = (f%, m ,j) = for all m = -1, 0, 1, 
the following estimate holds: 

(j,HMj)>C'\\jf. 

Similarly as in the discussion before Proposition [3J we note that / and f'Y l m are 
zero modes of Thus the proposition says that there are no other zero modes. 

Proof. We first note that is non-negative, since / is a minimizer of the one-polaron 
problem. We now argue that has exactly one negative eigenvalue. Indeed, since 
{j,Zr'j) = (j,H^j) if (f,j) = and since > 0, the variational principle implies 
that has at most one negative eigenvalue. On the other hand, the Euler-Lagrange 
equation for / implies that (/, iP-'f) = —A(f,Xff) < 0, which, again by the variational 
principle, means that L^ 1 ' has at least one negative eigenvalue. This proves the claim. 
We next recall Lenzmann's result [TJ, which states that 

kerL (1) = span{/'y lim : m = -1,0,1}. (3.3) 

We also note that commutes with angular momentum. Moreover, its essential 
spectrum starts at —fi/2 > and, therefore, (I3.3P implies that there is a constant 
d > such that 

(j,HWj) = (j,LWj)>d\\j\\* 

for all j with angular momentum I > 2 and for all j with angular momentum I = 1 
satisfying the additional constraint that (f'Yi tTn ,j) = for all m = —1, 0, 1. (Actually, 
the argument of [7J shows that the best constant d is achieved either for I = 1 or for 
1 = 2.) 

Since the function r appearing in the definition of is radial, the operator 
commutes with angular momentum. Therefore, the previous discussion reduces the 
proof of Proposition [5] to finding a lower bound on (j, H^'j) for radial j satisfying 
(f,j) = 0. In other words, we have to exclude the possibility that is a degenerate 
eigenvalue of restricted to I = 0. 

As an aside, before completing the proof, we show that there is a radial function R 
such that 

L {l) R = 

and 

(*,/) = \. 

This is also contained in [TJ, but we include the short proof for the convenience of the 
reader. We define fp{x) = (3 2 f((3x) and note that the Euler-Lagrange equation for / 
implies the following equation for fp, 



SYMMETRY OF BIPOLARON BOUND STATES — January 18, 2012 



11 



By differentiating this equation at ft — 1 we obtain L^^\p = xf^ = fif and thus 

R=4~o fp = 2f + x-Vf 

dp /3=1 

satisfies the desired equation and is obviously radial. The overlap (R, f) is computed 
by differentiating the identity J /^(x) 2 <ir = /3 at /3 = 1. 

Having found the function R we now conclude the proof by showing that is a 
simple eigenvalue of restricted to radial functions. Thus, let v be a radial function 
satisfying H^'v = and (/, v) = 0. We need to show that v — 0. By the expression 
for iff 1 ), 

L {1) v + (r,v)f = 0. 

We define v = v + /i _1 (r, v)R with R as constructed above. Thus L^v = 0. Since v is 
radial, Lenzmann's result (13. 3p implies that v = 0. From the orthogonality condition 
of v we infer that 

Thus, since (/, R) ^ 0, we have (r, v) = 0, and therefore v = v = 0, as claimed. This 
concludes the proof of Proposition [5l □ 

3.2. Proof of Proposition [4]. We consider the operator H in the subspace L 2 symrn {W & ) 
of all functions j G L 2 (M 6 ) satisfying j(x,y) = j(y,x). We decompose this space as 
^ mm (K 6 )= Ho ©H 1; where 



n = {af®f + f®g + g®f: a e C, {f,g) = 0} 



and 



Hi = spanj/}" 1 <g> spanj/j^ . 

Since hf = and since the rank one operators entering in the definition of H only 
involve terms of the form k £g> / + / £g> k (with k being either r or /), the operator H 
leaves both subspaces invariant and we can study it separately on each subspace. 

We observe that H coincides with L on Hi and that, moreover, the operator Xf^f 
vanishes on that space. Therefore, on Hi the operator L is just a sum of two one- 
body operators h. Since / is positive, it is the ground state of h and since the essential 
spectrum of h starts at —fi/2 > 0, h has a gap 5 > above zero. We conclude that 
L > 25 on Hi. 

We now turn to the space Hq. More precisely, we are only interested in the space 

Ho H ran P L = {/ <g> g + g ® f : (f,g) = {fY 1>m , g) = 0} . 

For j = f®g + g®f from this space we have 

(j,Lj)=2(g,(h-Ax f )g), 

and by Proposition [5j this is bounded from below by 2C"||g|| 2 = C"||j|| 2 . This com- 
pletes the proof of Proposition HI □ 
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Appendix A. Existence of a rotation invariant optimizer 

Our goal in this section is to prove Proposition [5J that is, the existence of a rotation 
invariant optimizer. The heart of the proof is the following lemma. 

Lemma 6. Assume thatip n 6 ff 1 (IR 6 ) is a sequence of normalized, rotation invariant 
functions such that y/p^ n has a weak limit ^fp in i/ 1 (IR 3 ). Then 



liminf £o[i/j n ] > e ( / pdx I . 



Proof. As in Step 2 in Subsection 12.11 we rewrite the energy in terms of the potential 
$n := pn * M^ 1 generated by p n := p^ n . By the Hoffmann-Ostenhof inequality [6] we 
obtain 

£t#n] = // (|V^ n | 2 + |V^n| 2 - $«(x)|Vn| 2 ~ <$> n(yM n?) dx dy 

J </R 3 xR 3 

+ -*- / |V$J 2 dx 



>( (|Vv^| 2 - $ n p n ) dx + ^- [ |V$ n | 2 cfe. (A.l) 

Of course, we may pass to a subsequence and assume that £0^71] has a finite limit. 
Thus we infer that the sequence $ n is bounded in if 1 (M 3 ) and, after passing to another 
subsequence if necessary, we may assume it has a weak limit $ in if 1 (R 3 ). Below we 
shall argue that 

lim / ($ nPn - $p) dx = . (A.2) 



n— >oo 



Once this is proved, we conclude from (lA.ip and the lower semi-continuity of the terms 
/ \^\fpn\ 2 dx and f\V$ n \ 2 dx that 

liminf £ [^ n ] > / (|Vv^| 2 -$p) da; + — / |V$| 2 dx. 

By the same argument as in (12. 3p and by the definition of the single-polaron energy 
we find that 

/ (IVVpI 2 -^) dx+^- [ \V<$>\ 2 dx> f \V^p\ 2 dx-D[p,p}>e( [ pdx] . 

The last inequality follows by scaling p. This is the claimed lower bound. 

We are left with proving (]A.2j) . First, by Sobolev embeddings $ n converges to $ 



weakly in L 6 and we have p E L 1 fl L 3 C L 6 / 5 . Thus, we only need to prove 

lim / $ n (p n - p) dx = . 
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At this point we use the spherical symmetry of p n and $„. Newton's theorem implies 
that < $ n (x) < 2|x| _1 . Thus, for any R > 0, 



/ ■ 

J\x\>R 



$n (Pn ~ P) dx 



which can be made arbitrarily small, uniformly in n, by choosing R large. On the 
other hand, 



$n (Pn - P) dx 



\x\<R 



< 



[ ^(VK+Vpf dx] 



1/2 



x\<R 



fa - y/pf dx 



1/2 



By Rellich-Kondrashov the second term on the right side tends to zero as n — > oo for 
every fixed R > 0. Moreover, the first term is bounded uniformly in n, since 

®l(yfa+Vp) 2 dx< ( [ & n dx) ([ ( v ^ + v fpfdx) . 

(The fact that the p n -term is uniformly bounded follows from the fact that p n is 
uniformly bounded in L 1 PI L 3 by Sobolev inequalities.) This concludes the proof of 
(IA.2j) . and therefore the lemma is proven. □ 

The following lemma is well known. We include the proof for the convenience of 
the reader. An instructive example to keep in mind is where ip n (x,y) = f{x)g n (y) + 
9n(x)f(y) where g n and / are if 1 (R 3 ) functions with disjoint support and where g n 
converges weakly to zero in H l {E?). 

Lemma 7. Assume thatipn is normalized and converges weakly to zero in if 1 (R 6 ) and 
assume that ^Jp^ n converges weakly in _f/' 1 (IR 3 ) to some v /p. Then f R3 pdx < 1. 



Proof. For any R > we write 



J\x\<R 



p^ n (x) dx 



\tp n (x,y)\ 2 dx dy + 



x\<R 



\ip n (x,y)\ 2 dxdy 



\y\<R 



= 1 - 
< 1 + 



\ip n (x, y)\ 2 dx dy + 



|x|>R,|y|>_R 



\ip n (x,y)\ 2 dx dy 



x\<R,\y\<R 



\ip n (x,y)\ 2 dxdy . 

'\x\<R,\y\<R 

Since if) n converges weakly to zero in if 1 (R 6 ), the Rellich-Kondrashov theorem implies 
that the last double integral on the right side tends to zero as n — > oo. Thus, again 
by the Rellich-Kondrashov theorem now applied to Jpip n , 



p(x) dx = lim 



p^ n {x)dx < 1 . 



'\x\<R "■- fuv J\x\<R 

Since this is true for any R > 0, we obtain the assertion. 



□ 
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Proof of Proposition GJ Let ip n be a minimizing sequence for e ^ mm with e s J mm < e . 
Since Eu[ip n ] is bounded, ipn is bounded in H l and, after passing to a subsequence, it 
converges weakly in H l (R 6 ) to a function ip. We decompose ip n = ip + ip n , where now 
ip n converges weakly to zero in /f 1 (lR 6 ), and note that 

i = ii^ir = ii^n 2 +ii^ir+o(i). (a.3) 

We claim that asymptotically, the energy of the minimizing sequence splits as follows, 

SuM = £v[1>, $»] + 8u$ n , $ n ] + o(l) . (A.4) 

Here $ n = p^ n * \x\ _1 , and we recall that the notation Eu[ip,<&] was introduced in 
( 12. 2p . The proof of ( 1A.4|) relies on the weak convergence in H 1 for the positive terms 
in the energy functional and on the fact that 

$nO)?AnO, y)ip{x, y) dx dy = 0(1) . 

The proof of the latter relation also uses the weak convergence in H 1 , together with 
the fact that $ n is bounded in L°°. The details are as in [U Eq. (5.13)] and are 
omitted. 

It follows from (TO) , together with (E3D and ( TA~3l that 

£um > iivir^tuvirv] + a - ihd^m + , (a.s) 

where v n = H^nll -1 ^- Here we interpret ^[H^H" 1 ^] as zero if ip = 0, and similarly 
for 8u[v n ]. 

Given (1A.5jl . it is easy to deduce that = 1. Indeed, we argue by contradiction 
and assume that ||-0|| < 1, which is the same, by (1A.3jl . as assuming that \\ip n \\ has 
a non-zero limit. Hence v n converges weakly to zero in H 1 (R 6 ). By the Hoffmann- 
Ostenhof inequality [6] the square roots of the corresponding densities a n = p Vn are 
bounded in if 1 (M 3 ) and hence, after passing to a subsequence, have a weak limit y/a 
in i^ 1 (M 3 ). Since all v n are rotation invariant, we learn from Lemmas [6] and [7] that 



lim inf Eu [v n ] > e ( / a dx) > e. 



On the other hand, we trivially have 1 ip] > e s J mm . (This also holds if ip = 

with our convention.) Thus (lA.5j) implies that 

e s r m = lim SvM > U\\ 2 e s u ymm + (1 - IH| 2 )e. 

n— >oo 

Since \\ip\\ < 1, this contradicts our assumption that e ^ mm < e . 

Thus we have shown that H^H = 1, and now ( IA.5I) implies that Eu[ip n ] > £u[ip]+o(l), 
from which we deduce that e ^J mm > Ejj[ip], that is, ip is a minimizer. This completes 
the proof of Proposition |2j □ 
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Appendix B. On the rotation invariant optimization problem 

In the next two propositions we present some interesting facts about the rotation 
invariant minimizers. This appendix is not needed in the rest of the paper. 

Proposition 8. IfU< 1, then e s J mm < e. On the other hand, if (12. 9ft has a rotation- 
invariant solution 0^6 iJ 1 (IR 6 ), then U < 4. In particular, e ^J mrn = e f or JJ > A. 

Proof. The first part of the proposition follows by a simple variational computation 
in the manner of Zhislin's theorem; see, e.g., [L2| Thm. 12.2]. We write ip(x,y) = 
[fR( x )VR(y) + ? 7-R( a; )/i?(l/)]/v / 2, where /# is the one-polaron function smoothly cut 
off at some large radius R, appropriately normalized (thereby making an error in the 
energy of the order o(l/R)) and where t)r is a normalized, radial function with support 
in the shell 10 R < \x\ < 11R. Thus the total energy is, by Newton's theorem 

e + o(l/R) + (U - 1) / \^)ldy< e + o(l/R) + . 

Jr3 \y\ 11R 

For sufficiently large R, this number is less then e. 

The second part is an adaptation of the N < 2Z + 1 theorem in [HI Sec. VI. a]. In 
the present situation N = 2 and the effective Z = J p v dx/U = 2/U. □ 

Our second result in this section concerns the correlation of the particles in the 
optimizing rotation invariant state. We state this in form of a rearrangement inequal- 
ity, where the rearrangement is defined as follows. We identify a rotation invariant 
function ifj(x,y) with a function u(r,s,t), where r = \x\, s = \y\ and t = x -y/\x\\y\. 
For fixed r,s > 0, we denote by u*(r, s, •) the unique, non- increasing function on 
[— 1, 1] which is equi-measurable with \u(r, s, -)| (in the sense of Lebesgue measure on 
[— 1, 1]). If u(r, s, t) came from a function ip(x,y), we also use the notation ip*(x,y) 
for u*(r,s,t). Then and ip* are equi-measurable and, in particular, 



Proposition 9. For any U > and any rotation invariant ip G i/ 1 (R 6 ) with 
one has £u\$*\ < £uW\ unless ip = ip* a.e. 

In particular, if ip is an optimizer for e s J mm , then for any fixed |x| and \y\, tp{x,y) 
is a non-increasing function oi x ■ y /\x\\y\. 

Proof. Since = p^* because of equi-measurability, the D[p, p] term in the energy 
functional does not change when ip is replaced by ip*. Moreover, the repulsion term 
improves (unless xjj = ip*) since, for fixed r, s > 0, 



(r 2 - 2rst + s 2 ) 1 / 2 J_ x (r 2 - 2rst + s 2 ) 1 / 2 

unless u(r,s,t) = u*(r,s,t). This is a simple rearrangement inequality proved as in 
[Til Thm. 3.4]. It uses the fact that (r 2 — 2rst + s 2 ) -1 / 2 is an increasing function of t. 
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For the kinetic energy we compute 



|V^| 2 

and similarly for V y ip, so that 

(iV.Vf + |V,Vf) dxdy 



du 



dr 



+ 



du 



dt 



8n 2 



oo poo /"l 



JO J-l 



du 



dr 



du 



ds 



'l-i 1 



du 



dt 



r 2 s 2 dt ds dr . 



We now argue that for any fixed r and s, one has 



and similarly for and 





du 


2 

dt > 


L 


du* 




dr 






dr 



dt. 





du 


2 pi 


du* 


/>-*>) 






dt 


dt> J (1-t 2 ) 


dt 



dt. 



Of course, these two inequalities will complete the proof of the proposition. For the 
proof of the first inequality we approximate §^(r, s, t) by h~ l (u{r + h, s, t) — w(r, s, t)). 
Then the t-integrals of the terms \u(r + h, s,t)\ 2 and \u(r, s,t)\ 2 do not change under 
symmetrization, whereas 



Re J u{r + h, s,t)u(r, s,t) dt < J 



u*{r + h, s, t)u*(r, s, t) dt 



by a simple rearrangement inequality. As h — > we obtain the first one of two claims. 
For the proof of the second inequality we consider the function g on § 2 , given in 
spherical coordinates (</?, 0) by g(ip, 6) = u(r, s, cos 9). Then 

2 



| V w g\ 2 du = 2n 





du 


/>-«■) 




dt 



dt. 



Since passing from u(r, s, •) to u*(r, s, ■) corresponds to the usual symmetrization of g 
on the sphere, the latter inequality follows from standard symmetrization results; see, 
e.g., []]. This finishes the proof of Proposition [9j □ 
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